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WHAT IS CLAIMED IS: 



1 . A method for performing convolution of an n -dimensional pattern 
p(x\ t x n ), with a smooth kernel d which is at least approximately equal to a 
separated-spline kernel, to generate data indicative of a convolution result r = 
Dp = d xp, where D is the circulant of d, and n is greater than one, said 
method including the steps of: 

(a) specifying the separated-spline kernel as k(x\, x n ) = k\(x\ )... k n (x n ), 
where kj has circulant K h k{x\, x n ) admits of an operator A = A\A 2 ...A m and Aj 
is an annihilation or flattening operator which operates on the circulant Kj of 
kernel kj in such a manner that AjKj is sparse (when Aj is an annihilation 
operator) or AjKj is almost everywhere a locally constant matrix (when Aj is a 
flattening operator); 

(b) processing pattern data indicative of the pattern p and kernel data 
indicative of the kernel k\ to generate additional data indicative of q\ = 
A\k\ x p for each row of the pattern p; 

(c) processing the additional data to backsolve A\r\ = q\ for said each 
row of the pattern to determine data indicative of n =k\ x p for said pattern, 
by performing a preliminary ignition step in which a small number of 
components of r\ are computed directly, and then determining all other ones of 
the components of r\ using a natural recurrence relation determined by the 
operator A\\ 

(d) processing the data indicative of r\ — k\ x p to generate data 
indicative of q 2 = A 2 k 2 x (k\ x p) T for each row of (k\ x p) T , where = 
(k\ x p) T denotes transposition of r\ ; and 

(e) processing the data indicative of q 2 = A 2 k 2 x r x T to backsolve A 2 r 2 = 
q 2 for said each row of (k\ x p) T thereby determining data indicative of r 2 = 
(k 2 x k\ x p) T for said pattern /?. 
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2. The method of claim 1, wherein n = 2, and step (e) includes the step 

of performing a transposition on the data indicative of r 2 = (fe x ki x p) to 
generate data indicative of r = (fe x fci x /?) for said pattern, where r is a close 
approximation to the convolution Dp . 

3. The method of claim 2, also including the steps of: 

(f) generating a pattern signal in response to convolved image data 
determined at least in part by the data indicative of r = (k 2 x k\ x p) for said 
pattern; and 

(g) causing the pattern signal to be incident on a set of optics, and 
projecting a pattern image from the optics onto a lithography target in response 
to the pattern signal. 

4. The method of claim 3, wherein the pattern signal is an optical signal. 

5. The method of claim 3, wherein the pattern signal is an electron beam 

signal. 

6. The method of claim 3, wherein the convolved image data employed 
in step (f) is said data indicative of r = (k 2 x k\ x p) for said pattern. 

7. The method of claim 1, wherein n is greater than 2, and also including 
the step of: 

(f) after determining the data indicative of r 2 T = (k 2 x k\ x pf for said 
pattern, performing;* iterations of steps (d) and (e), wherein during the z'-th 
iteration of steps (d) and (e), where i is a positive integer and/ = (n - 2), data 
indicative of <^ +2 = A l+2 k l+2 x # + i and of r,+2 T = (^+2 x k i+ \ x ... x ^ x ^) T 
are generated. 
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8. The method of claim 1, wherein the pattern data is stored in memory 
locations of a memory prior to step (b), and wherein step (d) includes the steps 
of: 

(i) storing the data indicative of r\ = k\ x p in the memory such that 
data indicative of each column of {k\ x p) occupies a row of the memory 
locations formerly occupied by the pattern data indicative of a row of p; and 

(ii) during generation of the data indicative of q 2 = A 2 k 2 x (k\ x p) T for 
said each row of {k\ x pf , reading data indicative of one said column of (k\ x 
p) from one said row of the memory locations. 

9. A method for performing convolution of a two-dimensional pattern 
p(x, y), with a smooth kernel d which is at least approximately equal to a 
separated-spline kernel, to generate data indicative of a convolution result r = 
Dp = dxp, where D is the circulant of d, said method including the steps of: 

(a) specifying the separated-spline kernel as k(x,y) = k\(x)k\(y) 9 where 



k l (x) = (R 2 -x 2 ) for \x\<R 




for \x\>R, 



and where the kernel k\(x) satisfies: 



k\(x) = k c (x) + k+(x) + k-(x) 9 



where k c (x) is a kernel 



k c (x) = (R 2 - x 2 ) for \x\<R 



= 0 



for x > 7?, 



k+(x) is a one-tailed Laplacian kernel 
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k + (x) = cd~^ for \x\>R 



= 0 for \x\<R 9 
and k.(x) is a one-tailed Laplacian kernel 



k_ (x) ^cd x for \x\ < -R 

= 0 for \x\ > R\ 

(b) processing pattern data indicative of the pattern p and kernel data 
indicative of the kernel k c , the kernel k+, and the kernel k. to generate 
additional data indicative of a cap convolution k c x p(x) for each row of the 
pattern p 9 a second convolution k+ x p(x) for said each row, and a third 
convolution L x p(x) for said each row, and adding together the additional data 
indicative of the cap convolution, the second convolution, and the third 
convolution for said each row to generate data indicative of q\ = k\ x p(x,y) = 
(k c x p(x,y)) + (k+ x p(x,y)) + (k. x p(x,y)) 9 wherein 

the data indicative of the cap convolution is generated by computing 
A c k c x p(x) for said each row, where A c is an annihilation operator defined as 

A c k c {x) := k c (x+3) - 3k c (x+2) + 2>k c (x+\) - k c (x) 9 and computing 
A c ' l (A c k c x p)(x) = k c x p(x) using the recursion relation 

k c (x+3) - A c k c (x) + k c (x) + 3£ c (x+2) - 3k c (x+l), 

the data indicative of the second convolution is generated by computing 
A+k+ x p(x), for said each row, where A+ is an annihilation operator defined as 

A+ k+(x) = k+(x + 1) - d x k+(x), and computing A+ l (A+k+ x p)(x) = k+ x 
p(x) using the recursion relation: 

(k+ x p)(x+l,y) = (A+k+ x p){x,y) + d~\k + x p)(x,y), 

and 
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the data indicative of the third convolution is generated by computing A. 
k. x p(x) 9 for said each row, where A. is an annihilation operator defined as 

A. k.{x) = k.(x - 1) - d x k.(x\ and computing A: l (AJc. x p)(x) = k.x p(x) 
by recursing from end to beginning using the recursion relation: 
5 (k. x p )(x - l,y) = (AJl x p)(x,y) + d\ k. x p)(x,y); 

(c) transposing the data indicative of (k\ x p(x,y)) to generate data 
indicative of (&i x p{x } y)) ; and 

(d) repeating step (b), this time by processing the data indicative of (k\ 
x p(x,y)) T , rather than the pattern data, with the kernel data, thereby generating 

10 data indicative of q\ = (k\ (x) x k\(y) x p(x,y)) . 

10. The method of claim 9, wherein the pattern data is stored in memory 
locations of a memory prior to step (b), and wherein step (c) is accomplished by 
storing the data indicative of (k } x p) in the memory such that data indicative 
15 of each column of (k\ x p) occupies a row of the memory locations formerly 

occupied by the pattern data indicative of a row of p(x,y). 



11. The method of claim 9, also including the step of: 
(e) transposing the data indicative of q\ T = (k\ x k\ x p) T to produce 
20 data indicative of r = k xp. 



12. The method of claim 11, also including the steps of: 

(f) generating a pattern signal in response to convolved image data 

determined at least in part by the data indicative of r = k x p; and 
25 (g) causing the pattern signal to be incident on a set of optics, and 

projecting a pattern image from the optics onto a lithography target in response 

to the pattern signal. 



13. The method of claim 12, wherein the pattern signal is an optical 

30 signal. 
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14. The method of claim 12, wherein the pattern signal is an electron 
beam signal. 

5 1 5. A method for performing convolution of a two-dimensional pattern 

p(x, y), with a smooth kernel d which is at least approximately equal to a 
separated kernel, to generate data indicative of a convolution result r(x T y) = Dp 
= dxp, where D is the circulant of d 9 said method including the steps of: 

(a) specifying the separated kernel as k(x,y)= k 2 {x)k 2 {y) J with k(x,y) 
1 0 approximating the kernel d, where 

k 2 (x) = 2cos 2 (7Dc/2R) for \x\<R, 

k 2 (x) = 0 for \x\>R; 

(b) processing pattern data indicative of each row of the pattern p and 
kernel data indicative of the kernel k 2 to 

(i) generate data indicative of Ak 2 x p for said each row, for -R 
15 -3 < x < n-l, where Ak 2 (x) x p{x) = B(p(x + R) + p(x + R + 1) -p(x 

-R)-p(x-R + 1)); and 

(ii) performing a recursion operation, using the recursion relation 
r(x + 2) =Ar(x) + B(r(x+l) — r(x)) + r(x - 1), and using the relations r(- 
R — 4) = r(-R — 3) = r(~R — 2) = 0 to ignite said recursion operation, to 

20 find ripe) = A' 1 (Ak 2 (x) x p(pc)) for said each row, thereby producing, at 

the end of a final iteration of step (ii), data indicative of 
k 2 x p; 

(c) transposing the data indicative of k 2 x p to produce data indicative of 
(k 2 x p) J ; and 

25 (d) repeating step (b), this time on rows of the data indicative of 

(k 2 x p) T rather than on the rows of the pattern data, thereby generating data 
indicative of convolution of all columns of k 2 x p with k 2 . 
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16. The method of claim 15, wherein steps (b), (c), and (d) are 
performed by a programmed processor, the processor is in an initial state at the 
start of step (b), and after step (d), the processor performs the additional steps 
5 of: 

(e) transposing the data indicative of convolution of all columns of 
k 2 x p with k 2 to produce data indicative of the convolution result r(x,y); and 

(f) after step (e), returning to the initial processor state. 

10 17. The method of claim 16, also including the steps of: 

(g) generating a pattern signal in response to convolved image data 
determined at least in part by the data indicative of the convolution result r(x,y); 
and 

(h) causing the pattern signal to be incident on a set of optics, and 

1 5 projecting a pattern image from the optics onto a lithography target in response 

to the pattern signal. 

18. The method of claim 17, wherein the pattern signal is an optical 

signal. 

20 

19. The method of claim 17, wherein the pattern signal is an electron 
beam signal. 

20. The method of claim 15, wherein the pattern data is stored in 
25 memory locations of a memory prior to step (b), and wherein step (c) is 

accomplished by storing the data indicative of k 2 x p in the memory such that 
data indicative of each column of (k 2 x p) occupies a row of the memory 
locations formerly occupied by the pattern data indicative of a row of p(x,y). 



30 



21. The method of claim 15, wherein B = 2sin 2 (7t/2L). 
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22. The method of claim 15, wherein 5=1. 

23. A computer system, comprising: 

a processor programmed with software for convolving first data 
indicative of an ^-dimensional pattern p(x h x n ), with a smooth kernel d at 
least approximately equal to a separated-spline kernel, to generate data 
indicative of a convolution result r = Dp = d x p, where D is the circulant of d, 
and n is greater than one, by processing second data, indicative of a separated- 
spline kernel k(x x , ...,*„)= k\(x\ )... k n (pc n ) which approximates d, where kj has 
circulant Kj t the kernel £(xi, x n ) admits of an operator^ =A\A 2 ..A n , where ,4,- 
is an annihilation or flattening operator which operates on the circulant Kj of 
kernel kj in such a manner that AjKj is sparse (when Aj is an annihilation 
operator) or AjKj is almost everywhere a locally constant matrix (when A 3 is a 
flattening operator), with said first data to generate third data indicative of q x = 
A\k\ x p for each row of the pattern p, and processing the third data to 
backsolve A x r x = q x for said each row of the pattern to determine fourth data 
indicative of r x = k\ x p for said pattern, and processing the fourth data to 
generate fifth data indicative of q 2 = A 2 k 2 x (k x x p) T for each row of (M x pf, 
where r\ T = (k\ x pf denotes transposition of r u and processing the fifth data 
to backsolve A 2 r T = q 2 for said each row of (k\ x p) T thereby determining sixth 
data indicative of r 2 = (k 2 x k\ x p) T for said pattern; 

at least one memory coupled to and accessible by the processor, and 
configured to store at least the first data, the second data, and the fifth data; and 

an input device coupled to the processor, wherein the processor executes 
the software in response to at least one control signal from the input device. 

24. The system of claim 23, wherein n=2, and the software causes the 
processor to perform a transposition on the sixth data to generate seventh data 
indicative of r =(k 2 x k\ x p) for said pattern, wherein the seventh data is said 
data indicative of the convolution result r = Dp = dx p. 
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25. A digital signal processor, including: 

digital signal processing circuitry configured to perform arithmetic 
operations on data including first data indicative of an ^-dimensional pattern 
p(x h x n ) , and second data indicative of a separated-spline kernel k(x x , 
x n )= h(x x )... k n (x n ) , where kj has circulant K 3 , the kernel k(x u *„) admits of 
an operator A =AiA 2 ..A„, where 4 is an annihilation or flattening operator 
which operates on the circulant Kj of kernel kj in such a manner that AjKj is 
sparse (when Aj is an annihilation operator) oxAjKj is almost everywhere a 
locally constant matrix (when A 3 is a flattening operator), and said kernel k is at 
least approximately equal to a smooth kernel d having circulant D, where n is 
greater than one; and 

a program control unit coupled to the digital signal processing circuitry, 
and configured to generate control bits in response to instructions, and to assert 
the control bits to the digital signal processing circuitry to cause the digital 
signal processing circuitry to perform convolution on the first data by 
processing the first data and the second data to generate third data indicative of 
qx = A x k x x p for each row of the pattern p, and processing the third data to 
backsolve A x r x = q x lox said each row of the pattern to determine fourth data 
indicative of r x = k x x p for said pattern, and processing the fourth data to 
generate fifth data indicative of q 2 - A 2 k 2 x (k x x pf for each row of (k x x p) 1 \ 
where r x T = (k x x p) 1 denotes transposition of r x , and processing the fifth data 
to backsolve A 2 r 2 - q for said each row of (k x x pf thereby determining sixth 
data indicative of r 2 = (k 2 x k x x pf for said pattern. 

26. The system of claim 25, wherein n = 2, and the control bits cause the 
digital processing circuitry to perform a transposition on the sixth data to 
generate seventh data indicative of r = (k 2 x k x x p) for said pattern. 

27. A lithography system, including: 

a digital signal processor, comprising digital signal processing circuitry 
and a program control unit coupled to the digital signal processing circuitry, the 
digital signal processing circuitry being configured to perform arithmetic 
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operations on data, including first data indicative of an ^-dimensional pattern 
p(x h x n ) , and second data indicative of a separated-spline kernel k(x h 
*n)= ki(x { )... k n {x n ) , where kj has circulant K h the kernel k(x h x n ) admits of 
an operator A = A x A 2 ...A n , where ^ is an annihilation or flattening operator 
which operates on the circulant K } of kernel kj in such a manner that AjKj is 
sparse (when A 3 is an annihilation operator) or AjKj is almost everywhere a 
locally constant matrix (when Aj is a flattening operator), and said spline kernel 
k is at least approximately equal to a smooth kernel d having circulant D, where 
n is greater than one, in response to control bits from the program control unit, 
to perform convolution on the first data by processing the first data and the 
second data to generate third data indicative of q x = A x k x x p for each row of 
the pattern/?, and processing the third data to backsolve A x r x = q x for said each 
row of the pattern to determine fourth data indicative of r x = k x x p for said 
pattern, and processing the fourth data to generate fifth data indicative of q 2 = 
A 2 k 2 x (k x x pf for each row of (k x x P ) T 9 where r x T = (k x x pf denotes 
transposition of r u and processing the fifth data to backsolve A 2 r 2 - q for said 
each row of (k\ xp) T thereby determining sixth data indicative of r 2 T = (k 2 x 
ki x pj 1 for said pattern, and performing a transposition on the sixth data to 
generate seventh data indicative of r 2 = (k 2 x k x x p) for said pattern; 

a pattern signal generator configured to generate a pattern signal in 
response to convolved image data determined at least in part by said seventh 
data; and 

a set of optics positioned so that the pattern signal is incident thereon, 
wherein the optics projects a pattern image on a lithography target in response 
to the pattern signal. 

28. The lithography system of claim 27, wherein the digital signal 
processor asserts the seventh data to the pattern signal generator, and the pattern 
signal generator is configured to generate the pattern signal in response to said 
seventh data. 
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29. The lithography system of claim 27, wherein the pattern signal is an 
optical signal. 

5 30. The lithography system of claim 29, wherein the optical signal is an 

optical beam and the set of optics is configured to focus the optical beam to 
generate a focused beam and to scan the focused beam relative to the 
lithography target. 



10 31. The lithography system of claim 27, wherein the pattern signal is an 

electron beam signal. 

32. The lithography system of claim 31, wherein the set of optics is 
configured to focus the electron beam signal to generate a focused electron 
15 beam and to scan the focused electron beam relative to the lithography target. 



33. An apparatus configured to convolve a first signal indicative of an n- 
dimensional pattern p(x\, x n ), with a second signal indicative of a separated- 
spline kernel k{x\, x n )= k\(x\ )... k n (x n ) which approximates a smooth kernel 

20 d, where kj has circulant K p the kernel k{x\, x n ) admits of an operator A = 

A\A 2 .,A n ^ where Aj is an annihilation or flattening operator which operates on 
the circulant K } of kernel kj in such a manner that AjK ; is sparse (when A 3 is an 
annihilation operator) or AjKj is almost everywhere a locally constant matrix 
(when Aj is a flattening operator), to generate a signal indicative of a 

25 convolution result r = Dp = d x p, where D is the circulant of d, and n is 

greater than one, by generating from the first signal and the second signal a 
third signal indicative of q\ = A\k\ x p for each row of the pattern /?, and 
processing the third signal to backsolve A\r\ = q\ for said each row of the 
pattern to generate a fourth signal indicative of r\ = k\ x p for said pattern, and 

30 generating from the fourth signal a fifth signal indicative of q 2 = A 2 k 2 x (k\ x 
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p) T for each row of (k x x pf, where r x J = (k x x p) T denotes transposition of r x , 
and processing the fifth signal to backsolve A 2 r 2 = qi for said each row of (k x x 
p) T thereby generating a sixth signal indicative of r 2 = (k 2 x k x x p) T for said 
pattern, whereby said sixth signal is at least approximately indicative of the 
convolution result r-Dp, 

34. The apparatus of claim 33, wherein the apparatus consists of 
electronic circuitry. 

35. The apparatus of claim 33, wherein the apparatus is a programmable 
gate array-based computing system. 

36. A computer-readable storage medium which stores instructions, 
wherein the instructions are executable by a computer, and wherein the 
instructions are such that the computer performs a method for convolving first 
data indicative of an 77 -dimensional pattern p(x\, x n ), with a smooth kernel d 
at least approximately equal to a separated-spline kernel k(x x , ...,*„)= k x (x x )... 
k n (x n ) 9 where kj has circulant K p the kernel k(x h x n ) admits of an operator^ = 
AiA 2 ..A n9 where A 3 is an annihilation or flattening operator which operates on 
the circulant Kj of kernel k 3 in such a manner that AjKj is sparse (when A } is an 
annihilation operator) or AjKj is almost everywhere a locally constant matrix 
(when Aj is a flattening operator), to generate data indicative of a convolution 
result r= Dp = d x p in response to executing said instructions, where D is the 
circulant of d, and n is greater than one, said method including the steps of: 

processing second data, indicative of the separated-spline kernel k(x h 
jc b )= Ati(jci )... k n (x n \ with said first data to generate third data indicative of q x = 
A\k\ x p for each row of the pattern p\ 

processing the third data to backsolve A x r x = q x for said each row of the 
pattern to determine fourth data indicative of r x = k x x p for said pattern; 
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processing the fourth data to generate fifth data indicative of q 2 = A 2 k 2 * 
(k x x p f for each row of (k } x p ) T ; where n T = x T denotes transposition 
of n; 

processing the fifth data to backsolve ^ 2 r 2 T = 92 for said each row of (k x 
x pf thereby determining sixth data indicative of r 2 T = (k 2 x k x x pf for said 
pattern; and 

performing a transposition on the sixth data to generate seventh data 
indicative of r =(k 2 x k x x p) for said pattern, whereby said seventh data is at 
least approximately indicative of the convolution result r = Dp. 



